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The bound states of fermions in the adjoint representation are of interest in supersymmetric models. We 
investigate the energy spectrum of the simplest - the gluino-gluon bound states - on several lattices in the 
quenched approximation. We use a static approximation for the gluino propagator. We find continuum limits of 
the splitting between the few lowest states, with the energy difference between the two lowest states of 354 ± 9 
MeV. 



Motivation 

The supersymmetric standard model enriches 
the conventional particle spectrum in all sec- 
tors. There is motivation to study states in 
strongly coupled N=l SQCD, where the parti- 
cle mass spectrum must be determined by non- 
perturbative techniques. As well their relevance 
for supersymmetric phenomenology, these bound 
states, in R-parity conserving theories, are of cos- 
mological interest as dark matter candidates. In 
addition to this, recent detector searches for such 
particles have sought to eradicate the window for 
theories with a light gaugino sector . 

In a low temperature regime we do not ex- 
pect gluino loops to contribute significantly to the 
propagation of a state. We study states contain- 
ing gluinos using quenched lattices and a static 
gluino propagator approximation. The system- 
atic error from quenching the quarks will be dom- 
inant. One of the simplest of these states is the 
bound states of gluons and a gluino. They can 
be constructed by measuring the propagation of 
a single heavy gluino coupled to the Yang Mills 
vacuum. 

The Static Approximation 

The contribution to the QCD action of a heavy 
fermion is given by 

The propagator for such fields is the Greens func- 
tion of the corresponding wave equation. In the 



limit Mq — > oo, we can discard space-like com- 
ponents of the covariant derivative and note that 
the group contribution is given by the phase evo- 
lution of the fermion as it moves. This provides 
us with a propagator for a static fermion; 

where P is the helicity projection operator and 
the sign depends on the direction of particle evo- 
lution. 

We see that this can be extended to encom- 
pass fermions in the adjoint representation sim- 
ply by demanding the replacement of the link 
variables Uf^ with elements of SU(3) in the ad- 
joint representation; the real 8x8 G^^ matrices. 
These we construct by observing that the combi- 
nation of fundamental elements with the gener- 
ators Tr(?7A''[/tA'') satisfy all the group require- 
ments. 

The Cubic Group 

We need to create a gluon field coupled to the 
static gluino propagator in order to measure the 
energy from the correlation C{t) with Euclidean 
time: 

< C S„ < X\t)\n >< n|A(0) > e"-^"*. 

Here X is the creation operator for the state, and 
n is a complete set of intermediate propagating 
states. 

For maximal signal and minimal contami- 
nation from other states we choose X to be 
irreducible representations of the cubic group 



0(8)^2 0^2(1 such that |< X\n >|2 is both 
large and orthogonal to other states. 

In the continuum limit states are labelled by 
J^'^' . In a lattice calculation the spatial rota- 
tional symmetry is broken down to that of O and 
we recover the spin content of O by subduction. 
The Spin content of O is given in Table . 

Table 1 

Irreducible representations of O 
Representation Dimension Spin Content 



Ai 1 0, 4, 6, 8... 

A2 1 3, 6, 7, 9... 

E 2 2, 4, 5, 6... 

Ti 3 1, 3, 4, 5, 5... 

T2 3 2, 3, 4, 5... 



We chose for simplicity the product of links 
in a square as the fundamental object in con- 
structing operators. This shape does not allow 
all J^^ combinations due to cancellations arising 
from the symmetries of the square. Combining 
the product of links to the adjoint propagator by 
a group generator, diagrammatically the correla- 
tion in a typical O representation looks like that 
in figure 0. 




Figure 1. Schematic gluelump sum 



Computational Details 

Measurement of objects containing static prop- 
agators are hampered by cumulative statistical 
errors from multiplying time links. We observe 
that adjoint links are even more sensitive to this 
effect. In order to make effective measurements at 
larger time to determine ground state masses, we 



employ a multi-hit technique. We run a heatbath 
algorithm locally on each time link and average 
over many samples to produce observables with 
significant variance reduction. 

We measured the correlations for a given state 
using four paths at both source and sink. These 
were constructed using two fuzzing levels on 
clover sums built from two sizes of squares. We 
then employed a variational technique the result- 
ing matrix to determine the lowest eigenvalue 
(See for example ) . We used a bootstrap analy- 
sis to determine the mass and statistical variation 
at a time separation of 2 or 3 lattice units, where 
by inspection of the error on the signal, a plateau 
had been reached. 

We measured correlations from all sites on var- 
ious quenched lattices, the statistics to date are 
shown in Table |^. The fuzzing levels and specific 
sizes were tuned according to the lattice spacing 
to give the best signal. 



Table 2 






Lattices used 


in calculation 






Size 


Number 


5.7 


8^ X 16 


20 


5.7 


12^ X 24 


20 


5.9 


123 X 24 


10 


6.0 


16^ X 48 


202 


6.2 


243 X 48 


60 



The Adjoint Static Spectrum and Contin- 
uum Limit Extrapolation 

Figure || shows the spectrum of states calcu- 
lated at /3 = 6.0. We have used a value ro/a{(3 = 
6.0) = 5.272 to scale the right hand axis inde- 
pendantly of the lattice spacing a. The points 
marked by circles are the ten measured represen- 
tations and are labelled as such. They are plotted 
assuming the lowest spin contained in the O rep- 
resentation. 

Note we confirm that the most symmetric state 
is not the ground state. 

If we consider a variational basis of large di- 
mension, and consider the restoration of rota- 
tional symmetry, then we may expect some of 
the higher eigenvalues of a variational analysis to 




1 2 3 

Lowest Spin Content, J 



Figure 2. Mass spectrum of states at f] — 6.0. 

correspond to higher spin content in the O rep- 
resentation. We plot some of the excited states: 
they are marked them by triangles and labelled 
according to their expected (continuum) spin con- 
tent. Where two states are expected to coincide in 
spin content they have been grouped. The spin 
assignments shown are in qualitative agreement 
with this degeneracy. With excited states added 
in this way we see clear Regge trajectories begin- 
ning to emerge in the spectrum. 

A direct extrapolation of the continuum mass 
spectrum is not feasible as the self energy of the 
static propogator has an ultraviolet divergence. 
The energy differences between states, however, 
have a well defined continuum limit and we ex- 
trapolate to a = the difference between the sec- 
ond, third and states and the lowest . 
These results are shown in table || and the ex- 
trapolation shown in figure The rp used to 
remove explicit a dependence from observables 
is also measured in the quenched approximation 
and introduces normalization errors 0(10%). 



Table 3 



Transition 


A{Moro) 


Mev 


Tf - - T+- 


0.898 ±0.022 


354 ±9 


T,- - T+- 


1.426 ±0.023 


562 ±9 


A++ - T+- 


2.667 ±0.068 


1053 ± 27 



As well as 0{a^) errors in measurements that 
we seek to remove in the continuum extrapola- 
tion, we also note that there may exist 0{'m?a?) 
corrections also. We see in figure that the 
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Figure 3. AMrp versus (a/rp)^, continuum limit 
extrapolation. 

and states are measured at different ener- 

gies. We expect, from their spin content, that 
in the continuum limit these two states become 
degenerate. The extrapolation to a = of these 
states is unclear with our present data, but a com- 
mon value is within statistical error. 

Conclusions 

The measured spectrum is in agreement with 
previous calculations and phenomenology sug- 
gesting that the 1+" and 1~~ states of the gluon 
field are the lowest lying. The lowest energy dif- 
ference is higher than a previous SU(2) calcula- 
tion which obtained a value of 200 ± 70 MeV [||. 

We note that this is a measurement of the J^'~' 
of the gluonic fields only. The lower lying states 
correspond to the simplest magnetic and electric 
modes of excitation of these fields. With explicit 
gluino spin included, these states split into degen- 



erate J: 



states. 



We see that finite size effects are under control 
and that extrapolation to the continuum limit of 
energy differences is well defined. 

The mass spectrum and structure of such states 
can act as phenomonological tool to aid under- 
standing of other gluonic states. 
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